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Abstract 



The Kelvin-Helmholtz instability in superfluids is discussed, based 
on the first experimental observation of such instability at the inter- 
face between superfluid 3 He-A and superfluid 3 He-B (R. Blaauwgeers, 
V. B. Eltsov, G. Eska et al, |gond-mat/0111343| ) [7]. We discuss why 
q ■ the Kelvin-Helmholtz criterion, the Landau critical velocity for nucle- 

ation of ripplons, and the free energy consideration all give different 
^ ■ values for the instability treshold. 

PACS: 47.20.Ma, 67.57.Np, 68.05.-n 



1. Classical Kelvin-Helmholtz (KH) instability. KH insta- 
bility belongs to a broad class of interfacial instabilities in liquids, 
gases, plasma, etc. It refers to the dynamic instability of the in- 
terface of the discontinuous flow, and may be defined as the instability 
of the vortex sheet. Many natural phenomena have been attributed 
to this instability. Most familiar of them are generation by wind of 
waves in the water, whose Helmholtz instability [Q] was first analyzed 
by Kelvin || , and flapping of sails and flags analyzed by Rayleigh Q 
(see recent experiments in [Q]). 
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Many of the leading ideas in the theory of instability were origi- 
nally inspired by considerations about inviscid flows. The corrugation 
instability of the interface between two ideal liquids sliding with along 
each other was first investigated by Lord Kelvin ||, gj. The critical 
relative velocity \v\ — V2\ for the onset of corrugation instability is 
given by 

1 PW2 -( Vl -v 2 f = ^F . (1) 



2 pi + P2 

Here a is surface tension of the interface between two liquids; p\ and 
P2 are their mass densities; and F is related to the external field 
stabilizing the position of the interface: typically it is the gravitational 
field 

F = g( P i - P2) ■ (2) 
The surface mode (ripplon) which is excited first has the wave vector 



'F/a , (3) 
and frequency 

, PWl +P2V2 ( ,s 

= k ■ • (4) 

Pi + P2 

The excited ripplon propagates along the interface with the phase and 
group velocity: u r i pp ion = {pivi + P2V%)/{pi + p%). 

However, among the ordinary liquids one cannot find an ideal one. 
That is why in ordinary liquids and gases it is not easy to correlate 
theory with experiment. In particular, this is because one cannot 
properly prepare the initial state - the planar vortex sheet is never in 
equilibrium in a viscous fluid: it is not the solution of the hydrody- 
namic equations if viscosity is finite. That is why it is not so apparent 
whether one can properly discuss its 'instability'. 

Superfluids are the only proper ideal objects where these ideas can 
be implemented without reservations, and where the criterion of in- 
stability does not contain viscosity. Recently the first experiment has 
been performed in superfluids, where the nondissipative initial state 
was well determined, and the well defined treshold has been reported 
[JU . The initial state is the nondissipative vortex sheet separating two 
sliding superfluids. One of the superfluids ( 3 He-A) performs the solid- 
body like rotation together with the vessel, while in the other one 
( 3 He-B) the superfluid component is in the so-called Landau state, 
i.e. it is vortex-free and thus is stationary in the inertial frame. The 



2 



threshold of the Kelvin-Helmholtz type instability has been marked 
by formation of vortices in the vortex-free stationary superfluid: this 
initially stationary superfluid starts to spin-up by the neighboring ro- 
tating superfluid. 

2. KH instability in superfluids at low T. The extension of 
the consideration of classical KH instability to superfluids adds some 
new physics. First of all, it is now the two- fluid hydrodynamics with 
superfluid and normal components which must be incorporated. Let 
us first consider the limit case of low T, where the fraction of the 
normal component is negligibly small, and thus the complication of 
the two-fluid hydrodynamics is avoided. In this case one may guess 
that the classical result (|I|) obtained for the ideal inviscid liquids is 
applicable for superfluids too, and the only difference is that the role 
of the gravity is played by the applied gradient of magnetic field H, 
which stabilizes the position of the interface between 3 He-A and 3 He-B 
in the experiment [?]: 



Here xa and xb are temperature dependent magnetic susceptibilities 
of the A and B phases. 

However, this is not the whole story. The instability will start 
earlier, if one takes into account that there is a preferred reference 
frame. It can be the frame of container, the frame of the crystal 
in superconductors, or even the frame where the inhomogeneity of 
magnetic field H is stationary. The energy of the excitations of the 
surface, ripplons, can become negative in this reference frame, and 
the surface modes will be excited, before the onset of the classical KH 
instability. 

Let us consider this phenomenon. We repeat the same derivation 
as in case of classical KH instability, assuming the same boundary 
conditions, but with one important modification: in the process of 
the dynamics of the interface one must add the friction force arising 
when the interface is moving with respect to the container wall. In the 
frame of the container, which coincides with the frame of the stable 
position of the interface, the friction force between the interface and 
container is 



F = (1/2)( X a(T)-xb(T))V(H 2 ) . 



(5) 



Ffriction 



(6) 
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where £(x, t) is perturbation of the position of the interface: 



z = zq + ((x, t) , ((x, t) = a s'm(kx — ut) . (7) 

We assume that the velocities v\ and v 2 are both along the axis x; 
the container walls are parallel to the (x, z)-plane; and the interface 
is parallel to the (x, y)-plane. 

The friction force in Eq.(^) violates the Galilean invariance in x- 
direction, which reflects the existence of the preferred reference frame 
- the frame of container. This symmetry breaking is the main reason 
of the essential modification of the KH instability. The parameter 
r in the friction force has been calculated for the case when the in- 
teraction between the interface and container is transferred by the 
normal component of the liquid due to Andreev scattering of ballistic 
quasiparticles by the interface ||. The friction modifies the classical 
spectrum of surface modes: 

(uj \ 2 fu \ 2 F + k 2 a _u 

Pl {k- Vl ) +p2 {k- v V =^—- iT k' (8) 

or 

k Pi + P2 

, 1 lF + k 2 a ptp 2 " ~ f . 

±-==J iT- (v 1 -v 2 y, 9 

VPl + P2 V « « Pi + P2 

where v\ and v 2 are the velocities of superfluid components of the 
liquids with respect to the container walls. 

For r = the spectrum of ripplons acquires the imaginary part, 
Imw(fc) ^ 0, at the classical treshold value in Eqs.((l|) and (||). How- 
ever, the frame-fixing parameter V leads to essentially different result: 
The imaginary part of frequency becomes positive Ymuj{k) > first 
for ripplons with the same value of the wave vector, as in Eq.(||), 
but the ripplon frequency is now u = and its group velocity is 
v group = dco/dk = 0. The critical ripplon is stationary in the reference 
frame of the container, as a result the onset of instability is given by 

\pivj + ~piv\ = V^F . (10) 

This criterion does not depend on relative velocities of superfluids, but 
is determined by velocities of each of the two superfluids with respect 
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to the container (or to the remnant normal component). Thus the 
instability can occur even if two liquids have equal densities, p\ = p2, 
and move with the same velocity, v\ = V2- This situation is very simi- 
lar to the phenomenon of flapping flag in wind, discussed by Rayleigh 
in terms of the KH instability - the instability of the passive de- 
formable membrane between two distinct parallel streams having the 
same density and the same velocity (see latest experiments in Ref . j| ) . 
In our case the role of the flag is played by the interface, while the 
role of the flagpole which pins the flag (and thus breaks the Galilean 
invariance) is played by the container wall. 

Note that in the limit of the vanishing pinning parameter r —* 
the Eq.([l(]) does not coincide with the classical equation ([[]) obtained 
when there is no pinning, i.e. when T is exactly zero. Such difference 
between the limit and exact cases is known in many area of physics. 
In classical hydrodynamics the normal mode of inviscid theory may 
not be the limit of a normal mode of viscous theory ||. Below we 
discuss this difference for the case of KH instability in superfluids. 

3. KH instability and modified Landau criterion. Let us 
first compare both results, with no pinning (r = 0) and for vanishing 
pinning (r — ► 0), with the Landau criterion. According to Landau, a 
quasiparticle is created by the moving superfluid if its velocity with 
respect to the container wall (or with respect to the body moving in 
superfluid) exceeds 



Let us recall that the energy E(p) here is the quasiparticle energy 
in the reference frame moving with the superfluid vacuum. In our 
case there are two superfluids moving with different velocities. That 
is why there is no unique superfluid comoving frame, where E(p) can 
be uniquely determined. Such frame appears only in particular cases, 
when either v% = V2, or if instead of the interface one considers the 
free surface of a single liquid (i.e. if p2 = 0). In these particular cases 
the Landau criterion in the form of Eq.(|Tl|) must work. The energy 
spectrum of the ripplons on the interface between two stationary fluids 
(or on the surface of the single liquid) is given by Eq.(|9|) with v\ = 



E{p) 




^Landau = mill 



P 



v 2 = T = 0: 



L0 2 (k) 



1 F + k 2 a 
Pi+ P2 k 
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This gives the following Landau critical velocity: 

L0 2 (k) 2 r— 
^Landau = mm ~^2~ = T ^ F<J . (13) 

k Pi + P2 

This coincides with the Eq.(|To|) if v% = V2, or if p2 = 0. But this does 
not coincide with the classical KH result: the latter is obtained at 
T = when the interaction with the reference frame of the container 
is lost, and thus the Landau criterion is not applicable. 

In the general case, when neither of the two conditions (vi = 1)2, 
or p2 = 0) fulfils, the Landau criterion must be reformulated: the in- 
stability occurs, when the frequency of the surface mode in the frame 
of the container crosses zero for the first time: u(k;vi,V2) = 0. In- 
spection of Eq.Q with T = shows that for k = ko the spectrum 
with negative square root touches zero just when the treshold (|l0| ) is 
reached. Thus the Landau criterion in its general formulation coin- 
cides with the criterion of instability obtained for the case of nonzero 
friction force. As distinct from the Landau criterion in the form of 
( |Tl"D valid for a single superfluid velocity, where it is enough to know 
the ripplon spectrum in the frame where the superfluid (s) is (are) 
at rest, in the general case one must calculate the ripplon spectrum 
u(k;vi,V2) for the relatively moving superfluids. 

4. Matching zero-pinning and vanishing-pinning regimes. 
The difference in the result for onset of KH instability in the two 
regimes - with T = and with r / - disappears only in the case 
when two superfluids move in such a way that in the reference frame 
of container the combination p±vi + P2V2 = 0. In this arrangement, 
according to Eq.(||), the frequency of the ripplon created by classical 
KH instability is zero in the container frame. Thus at this special 
condition the two criteria, zero pinning (|l|) and vanishing pinning ([!(]) , 
must coincide; and they really do. 

If pi^i + P2V2 7^ 0, the crossover between the zero pinning regime 
and the regime of small pinning occurs by varying the observation 
time. Let us consider this on the example of the experimental set-up 
0] with the vortex-free B- phase and the vortex-full A-phase in the 
rotating vessel: In the container frame one has vi = v s a = , V2 = 
v s # = —12 x r; the densities of two liquids, 3 He-A and 3 He-B, are 
the same with high accuracy: Pa = PB = P- In the non-zero pinning 
regime the instabilty occurs at the boundary of the vessel, where the 
velocity of the 3 He-B is maximal, when this maximal velocity reaches 
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the value: 

V 2 C = -Vfr = -4h = 2^Landau ■ (14) 
P 2 

This velocity is by y/2 smaller than that given by classical KH equation 
(|l]) for the zero-pinning regime. On the other hand it is by y/2 larger 
than the Landau criterion in the form of Eq. ([Til ) j but coincides with 
Landau criterion properly formulated for two superfluids. 

From Eq.(|8|) it follows that slightly above this treshold the incre- 
ment of the exponential growth of the interface perturbation is 

Im w(feo) = ( — - l) , at v sB -»c«" c . (15) 

In the vanishing pinning limit r — > the increment becomes small and 
the discussed instability of the surface has no time for development 
if the observation time is short enough. It will start only at higher 
velocity of rotation when the classical treshold of KH instability, v kh 
in Eq. (]l]), is reached. Thus, experimental results in this limit would 
depend on the observation time - the time one waits for the interface to 
be coupled to the laboratory frame and for the instability to develop. 
For sufficiently short time one will measure the classical KH criterion 
(|]), while for the sufficiently long observation time the modified KH 
criterion (^) will be observed. 

5. Thermodynamic instability. Let us now consider the case of 
nonzero T, when each of the two liquids contain superfluid and normal 
components. In this case the analysis requires the 2 x 2-fluid hydrody- 
namics. This appears to be rather complicated problem, taking into 
account that in some cases the additional degrees of freedom related to 
the interface itself must be also added. The two- fluid hydrodynamics 
has been used for investigation of the instability of the free surface 
of superfluid 4 He by the relative motion of the normal component of 



the liquid with respect to the superfluid one [10]. We avoid all these 
complications assuming that the viscosity of the normal components 
of both liquids is high, as it actually happens in superfluid 3 He. In 
this high-viscosity limit we can neglect the dynamics of the normal 
components, which is thus fixed by the container walls. Then the 
problem is reduced to the problem of the thermodynamic instability 
of the super flow in the presence of the interface. 

We start with the following initial non-dissipative state correspond- 
ing to the thermal equilibrium in the presence of the interface and 
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superflows. In thermal equilibrium the normal component must be 
at rest in the container frame, v n x = v„2 = 0, while the supernuids 
can move along the interface with velocities v s i and v S 2 (here the 
velocities are in the frame of the container). 

The onset of instability can be found from free-energy consider- 
ation: When the free energy of static perturbations of the interface 
becomes negative in the frame of the container, the initial state be- 
comes thermodynamically unstable. The free-energy functional for 
the perturbations of the interface in the reference frame of the con- 
tainer is determined by 'gravity', surface tension, and perturbations 
v sl = V<3?i and v S 2 = V<1>2 of the velocity field caused by deformation 
of the interface: 

HQ = \J dx(Fe + °{d x Cf + 

+ J dzp slik v % sl v k sl + dzp s2 ikV l s2 v k s ^j ■ ( 16 ) 

For generality we discuss anisotropic supernuids, whose supernuid 
densities are tensors (this occurs in 3 He-A). The velocity perturbation 
fields v s k = V<£fc) obeying the continuity equations di(p lk v s k) = 0, 
have the following form: 

$i(x,z < 0) = Aie klZ coskx, , . 

<&2(x,z > 0) = A 2 e~ k2Z coskx, ^ ' 

Pslzki = Pslxk 2 , Ps2zk\ = Ps2xk 2 . (18) 

The connection between the deformation of the surface, C( x ) = a s ^ n ^ x t 
and the velocity perturbations follow from the boundary conditions. 

Because of large viscosity of the normal component it is clamped 
by the boundaries of the vessel. Then from the requirement that the 
mass and the heat currents are conserved across the wall, one obtains 
that the supernuid velocity in the direction normal to the wall must 
be zero: v s i • n = v S 2 • n = 0. This gives the following boundary 
conditions for perturbations: 

d z ^ = v sl d x (, d z ®2 = v s2 d x ( . (19) 



Substituting this to the free-energy functional (16), one obtains the 



quadratic form of the free energy of the surface modes 
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X \ F + k 2 a - k (y/PsxlPszlvli + y/ Psx2Psz2V 2 s2j) (20) 

This energy becomes negative for the first time for the mode with 
k = {F/a) 1 / 2 when 

2 (\/PsxlPszlV 2 sl + ^ Psx2Psz2V 2 s2 ) = VcrF • (21) 

This is the criterion ( JTo| ) for the non- zero pinning regime extended 
to finite temperatures. Eq.(^l|) transforms to Eq.(|l~0|) when T — > 0: 
The normal components of the liquids disappear and one has p sx \ = 
Pszi = Pi and p sx2 = Psz2 = P2- 

6. Nonlinear stage of instability. Eq.(21) is in excelent agree- 
ment with the onset of the surface instability measured in experiment 
[|7j. The onset of instability is marked by the appearance of the vortex 
lines in 3 He- B which are monitored in NMR measurements. This 
demonstrates that vortices appear in the nonlinear stage of this KH 
instability. 

The precise mechanism of the vortex formation is not yet known. 
One may guess that the A-phase vorticity is pushed by the Magnus 
force towards the vortex-free B-phase region jy]]. When the potential 
well for vortices is formed by the corrugation of the interface (see 
Figure) , 

the vortices are pushed there and enhance further the growth of 
the potential well, until it forms the droplet of the A-phase filled 
by vorticity. The vortex-full droplet propagates to the bulk B-phase 
where it relaxes to the singular vortex lines of 3 He-B. 

Under the conditions of the experiment the nucleation of vortices 
leads to decrease of the B-phase velocity below the instability treshold, 
and the vortex formation is stopped. That is why one may expect that 
the vortex-full droplet is nucleated during the development of the in- 
stability from a single seed. The size of the seed is about one-half of the 
wavelength Ao = 2-Tr/ko of the perturbation. The number of the cre- 
ated vortices is found from the circulation of superfluid velocity carried 
by the piece of the vortex sheet of size Ao/2, which is determined by 
the jump of superfluid velocity across the sheet: k = \v s B — v s a| Ao/2. 
Dividing this by the circulation quantum kq of the created B-phase 
vortices one obtains the number of vortices produced as the result of 
the growth of one segment of the perturbation: 

N = J L ^v<+o (22) 
K 2ko 
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Figure 1: Possible scenario of vortex formation by Kelvin-Helmholtz insta- 
bility of the AB interface. 
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It is about 10 vortices per event under condition of the experiment, 
which is in a good agreement with the measured number of vortices 
created per event [f?J. This is in favour of the droplet mechanism of 
vortex formation. 

Probably, the experiments on KH instability in superluids will al- 
low to solve the similar problem of the non-linear stage of instability 
in ordinary liquids (see, for example, Ref. 

The vortex formation by surface instability is rather generic phe- 
nomenon. This mechanism has been discussed for vortex formation 
in the laser manipulated Bose gases |l3|, [l4j]. It can be applicable to 
different kinds of interfaces, and under very different physical condi- 
tions. In particular, vortices can be generated at the second order 
phase boundary between the normal and the superfluid phases 
Such an interface naturally appears at the rapid phase transition into 
the superfluid state [16|. The instability of the free surface of super- 
fluid under the relative flow of the normal and superfluid components 
of the same liquid has been recently reexamined by Korshunov [[17|]. 
He also obtained two criteria of instability: for zero and nonzero values 
of the viscosity of the normal component of the liquid. 

I thank R. Blaauwgeers, V. B.Eltsov, N. Inogamov, N. B.Kopnin, 
S. E. Korshunov, M. Krusius, E. A. Kuznetsov, and E.V. Thuneberg for 
fruitful discussions. This work was supported by ESF COSLAB Pro- 
gramme and by the Russian Foundations for Fundamental Research. 
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